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ABSTRACT 

GOODS 850-5 is a hyperluminous radio-faint submillimeter source in the GOODS-N. Although it is gen- 
erally agreed that GOODS 850-5 is at a high redshift z > 4, its exact redshift is unknown. While its stellar 
SED suggests z ~ 6, its radio/FIR SED suggests a lower redshift of z ~ 4. To better constrain its stellar SED 
and redshift, we carried out nano-Jansky sensitivity ultradeep NIR observations between 1.2 and 2.1 /im with 
the HST and the 8 m Subaru Telescope. Even with such great depths we did not detect GOODS 850-5, and 
the results show that it is an extremely curious source. Between the K s and 3.6 /im bands its spectral slope 
is > 3 x that of an ERO, and the flux ratio between the two bands is > 8 x that of Lyman breaks. It is quite 
challenging to explain this unusually red color without a Lyman break (which would imply z > 17). It requires 
a large amount (M* ~ 10 1L5 M Q ) of reddened old stars at z ~ 6, coexisting with an even more extinguished vi- 
olent ~ 2400^1-400 Mq yr -1 starburst, which does not have any associated detectable rest-frame UV radiation. 
We discuss the discrepancy between the NIR and radio/FIR photometric redshifts. We conclude that GOODS 
850-5 is at least at z > 4 and is more likely at z > 6. We describe the unusual properties of GOODS 850-5, 
including its SED and formation history, and we discuss the implications of such massive z > 6 galaxies. 
Subject headings: cosmology: observations — galaxies: evolution — galaxies: formation — galaxies : star- 
burst — infrared: galaxies — submillimeter 



1. INTRODUCTION 

The Rayleigh-Jeans portion of the dust spectral energy dis- 
tribution (SED) of IR-luminous galaxies produces a strong 
negative ^-correction and makes the observed submillime- 
ter flux of such galaxie s almost invariant at z > 1 to z ~ 10 
(iBlain & Longairll993l) . Thus if there are IR-luminous galax- 
ies at high redshift, observations at submillimeter wavelengths 
are a powerful way to find them. However, to date, all but one 
of the identified submillimeter galaxies (SMGs), other than 
those around luminous quasars, are at redshifts lower than 4, 
likely because of the limited resolution of the current submil- 
limeter instruments and the limited sensitivity of the current 
radio instruments, which are needed to locate the sources. 

The Submillimeter Common-User Bolometer Array 
(SCUBA) on the single-dish James Clerk Maxwell 
Telescope (JCMT) resolved 20%-30% of the sub- 
millimeter extragalactic background light (EBL) 
into point sources brighte r than ~ 2 mJy at 850 
um (ISmail Ivison & Blainl 11997b IBarger et al l 1 1998 : 
Hughes et all 119981; iBarger. Cowie. & Sanders! Il999 : 
Eale s et al. I~ ll999l) . Because of the low resolution of 
JCMT (~ 15" at 850 /im), identifications of the submil- 
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Telescope, obtained at the Space Telescope Science Institute, which is op- 
erated by the Association of Universities for Research in Astronomy, Inc., 
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2 Based in part on data collected at Subaru Telescope, which is operated 
by the National Astronomical Observatory of Japan. 

3 Jansky Fellow 
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NM 87801. The NRAO is a facility of the National Science Foundation op- 
erated under cooperative agreement by Associated Universities, Inc. 
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Charter Street, Madison, WI 53706 

6 Department of Physics and Astronomy, University of Hawaii, 2505 Cor- 
rea Road, Honolulu, HI 96822 

7 Institute for Astronomy, University of Hawaii, 2680 Woodlawn Drive, 
Honolulu, HI 96822 



limeter sources have to assume the radi o-far-infrared ( FIR) 
correlation in local galaxies (see, e.g.. ICondonlll992l) and 
rely on radio interferometry to pinpoint the location of 
the submillimeter emission. Optical spectroscopy of radio 
identified SMGs shows that they are ultraluminous sources 
(> 10 12 L©, corresponding to a star formation rate of 1 2 — 1 3 
Mq yr" 1 ) peaking at z ~ 2-3 and that the y dominate the 
total star f ormation in this redshift range (Chapm an et al.l 
l2003ll2005l) . However, the positive ^-correction of the radio 
synchrotron emission makes the radio wavelength insensitive 
to high-redshift galaxies, and radio observations can only 
identify 60%-70% of the blank-field submillimeter sources 
(IBarger. Cowie. & Richards! l2000t llvison et al l 120021) . The 
radio unidentified SMGs are commonly thought to be at 
redshifts higher than the radio detection limit (typically 
z > 3^1) but there has been a lack of direct evidence for such 
a high-redshift radio-faint submillimeter population. To date, 
there is only one s pectroscopically confirmed SMG at z > 4 
dCapak et all 120081, z = 4.547). Luminous radio-faint SMGs 
contribute ~ 10% to the submillimete r EBL measured b y 
the COBE satellite dPuget et al.l 119961: iFixsen et al.1 119981) . 
If most of these galaxies are indeed at high redshifts, then 
this implies a large amount of star formation in massive 
high-redshift systems. 

With recent developments in submillimeter interferometry, 
it is now possible to directly locate submilli meter sources 
withou t relying on radio interferometers. Young er" et al.1 
(120071) observed a sample of SMGs with the Submillime- 
ter Array (SMA). Several of their SMA detections are radio 
faint and consistent with being at redshifts higher than the 
radio identified SMGS. We have also been carrying out a pro- 
gram specifically targeting rad i o-fain t submillimeter sources 
with the SMA. In IWang et all d2007l hereafter W07) we re- 
ported our first identification in this program, GOODS 850-5. 
GOODS 850-5 was detected in our JCMT SCUBA jiggle-map 
survey of the Great Observatories Or igins Deep Survey-North 
(GOODS-N, Giavali sco et~ai1l2004l) with an 85 fjm flux of 
12.9 ±2.1 mJy (Wa ng. Cowie. & Bargerll2004l) . It was also 
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FIG. 1 . — Our new ultradeep NIR J, F160W, and K s band images of GOODS 850-5. Each panel has a 24" width. North is up. The three images have identical 
surface brightness scales (in /„). The 2" diameter circles mark the SMA position from W07, which has an uncertainty of Of'2. 



detected in the c ombined jiggle a nd scan map of the GOODS- 
N (GN 10, see iPope et al. 2006 and references therein). It is 
the second brightest submillimeter source in our jiggle-map 
catalog of the GOODS-N and has a total IR luminosity of 
~ 2 x 10 13 Lq. It did not have a 5 a radio counter part in the 
deep V ery Large Array ( VLA) 1 .4 GHz catalogs of lRichardsl 
(2000) and Big gs & Ivisonl d2006l) . 

In W07, an SMA 880 pm detection of GOODS 
850-5 was obtained and its counterpart was found to 
be extremely faint in the optical a nd near-IR (NIR). 
iDannerbauer. Walter. & Morri son (2008, hereafter D08) soon 
followed up with new 1.25 mm and 20 cm detections of 
GOODS 850-5. It turns out that t he correct Spitzer counter- 
part had already been suggested by Pop e et all ((2006) but this 
identification was not confirmed and the high redshift of the 
source was not realized until the accurate position was ob- 
tained by W07. However, the exact redshift of GOODS 850-5 
is unclear. Both W07 and D08 found that the NIR SED is con- 
sistent with a galaxy at z ~ 6, and the submillimeter and radio 
SED is consistent with z ~ 4. While both groups agree that 
this is a high-redshift SMG, W07 favors z ~ 6 but D08 favors 
z ~ 4 based on the radio-FIR correlation. The z ~ 6 redshift 
suggested by W07 is based on a featureless power-law con- 
tinuum in the IRAC bands and a non-detection in a relatively 
shallow K s band image. However, it is fair to say that this pho- 
tometric redshift is not a secure one, which would require the 
detection of at least one prominent spectral feature. To obtain 
a better constraint on the redshift, we carried out ultradeep 
NIR imaging in the J, H, and K s bands, hoping to detect the 
stellar continuum at < 2 fim. However, we found instead that 
the SED of GOODS 850-5 at 1-3 ^m is quite weird — it is not 
detected even at a 5 nJy sensitivity at 1 .6 /im (compared with 
its ~ 1 /iJy flux at 3.6 fim). Explaining this is challenging, 
but our analyses show that it greatly strengthens the previous 
W07 suggestion of z > 6. 

In this paper we present our new ultradeep NIR imaging 
observations of GOODS 850-5 and complete analyses of the 
likelihood of its redshift. The new observations and data used 
in this paper are described in § [2] Photometric redshift anal- 
yses in the optical, NIR, FIR, and radio are described in § [3] 
The mass and age of the stellar population and the star forma- 
tion rate of GOODS 850-5 are estimated, and the nature of this 
galaxy is discussed in §H The implications for galaxy evolu- 
tion are discussed in §|5] We summarize and make some final 
remarks in § [6] The cosmological parameters adopted in this 



paper are Ho = 71 km s 1 Mpc 51a = 0.73, and VIm = 0.21 . 

2. NEW NIR IMAGING AND EXISTING DATA 

2.1. HST NICMOS F160W Imaging 

We observed GOODS 850-5 with the Near Infrared Cam- 
era and Multi-Object Spectrometer (NICMOS) on the Hubble 
Space Telescope (HST) in Cycle 16. We used the NIC3 cam- 
era and the F160W filter to obtain the highest sensitivity. The 
observations were made in November and December 2007, 
containing a total of 16 HST orbits in four visits. In each orbit, 
we made only one exposure that is as long as possible (~ 47 
minutes) in order to minimize the contribution of read noise. 
The exposures were dithered to provide 0.5 pixel sampling. 

The data reduction was carried out in the Interactive Data 
Language (IDL) environment. We started the reduction with 
the standard pipeline-calibrated images provided by the HST 
Archive, which are flux calibrated and have instrumental sin- 
gatures removed. We applied a background subtraction to 
each exposure by fitting the object masked image with a 
smooth polynomial surface. The brightest cosmic ray hits 
and hot pixels were removed in each image with a spatial 
sigma filtering. The offsets between images were determined 
by measuring the positions of the detected objects in the cen- 
tral part of the field of view. We did not attempt to correct for 
the optical distortion of NIC3. The 16 images were then driz- 
zled to a common grid of 0."1 per pixel (two times finer than 
the original pixel of NIC3) and combined to form a deep im- 
age. The abs olute astrometry was m atched to the GOODS-N 
ACS catalog dGiavalisco et alj2 004) but corrected f or the "4 
offset in the ACS frame to match the radio frame (iRichardsl 
2000). Before the images were combined, another sigma 
filtering was applied to drizzled pixels in the same coordi- 
nate grids to remove fainter cosmic ray hits. When the im- 
ages were combined, each pixel was weighted by its exposure 
time and each image was inverse-weighted by the background 
brightness to achieve the highest S/N. The reduced image is 
shown in Figure Q] The final image has an extremely deep 1 
a sensitivity of 4.9 nJy. Surprisingly, even with such a great 
depth, we are not able to detect GOODS 850-5. 

2.2. Subaru MOIRCS J and K s Imaging 

We obtained extremely deep J and K s images of GOODS 
850-5 with the Mul ti-Object InfraRed C amera and Spectro- 
graph (MOIRCS, llchikawaet ail 120061) on the 8 m Subaru 
Telescope on Mauna Kea. MOIRCS contains two 2k x 2k 
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FIG. 2. — Multiwavelength images of GOODS 850-5. Each panel has a 24" width. North is up. The SMA position is labeled with 2" diameter circles and 
has an uncertainty of 0"2. The color codes in the optical, NIR, and IRAC images are labeled in the pictures. Grayscale images have negative color scales. The 
MOIRCS and NICMOS images are from this work and the VLA image is from D08 and G. Morrison (in preparation). The rest are from W07. 



HAWAII2 detectors, covering a field of view of ~ 4' x 7' with 
a pixel scale of 0." 117. The images used here were taken in 
two large MOIRCS imaging campaigns in the GOODS-N led 
by Japanese groups and our group based in Hawaii between 
December 2005 and January 2008. All the J band exposures 
and approximately half of the K s band exposures were made 
by Japanese investigators and were obtained from the Sub- 
ar u public archive . Parts of the Japanese data were published 
in lKajisawa et al.1 d2006). The rest of the K s exposures were 
made by our group and will be published elsewhere. The ma- 
jority of the observations were performed under photometric 
conditions with excellent seeing of 0"25-0"6. A very small 
fraction of the data has a large extinction of > 0.5 magnitude 
or poor seeing of > 0."7, and they were excluded in this work. 

The MOIRCS images were reduced with the IDL based 
SIMPLE Imaging and Mosaicking Pipeline (SIMPLE, W.-H. 
Wang 2008, in preparation 8 ). Images within a dither set (typ- 
ically 20-30 minutes in length) were flattened using an iter- 
ative median sky flat in which a simple median sky was first 
derived to flat the images and then a second median sky was 
derived by masking all the detected objects using the flattened 
images. After the images were flattened, the residual sky 
background was subtracted with a smooth polynomial sur- 
face. The brightest cosmic ray hits were removed by a spatial 
sigma filtering in each flattened image. MOIRCS produces 
almost nearly circular fringes in roughly half of the images. 
The fringes were modeled in polar coordinates where they are 
nearly perfect straight lines and were subtracted from the im- 



ages in the original Cartesi an coordinates. 

The package SExtractor dBertin & Arnoutsll 19961) was used 
to measure object positions and fluxes in each flattened, sky 
subtracted, and fringe removed image in a dither set. The 
first-order derivative of the optical distortion function was de- 
rived by measuring the offsets of each object in the dither se- 
quence as a function of location in the images. Absolute as- 
trometry was obtained by matching the detected objects to a 
reference catalog constructed with brighter and compact ob- 
jects in the GOODS-N ACS catalog (matched t o the radio 
frame) and the GOODS-N SuprimeCam catalog ( Cap ak et al.l 
2004). The images were then warped directly from the raw 
frames to a common tangential sky plane with a sub-pixel 
accuracy. All projected images were weighted by their sky 
transparencies, exposure times, inverse background bright- 
nesses, and pixel-to-pixel efficiencies (i.e., flat field) to ob- 
tain optimal S/N. The weighted images were then combined 
to form a large mosaic. When images from a dither set were 
combined, a sigma filter was applied to pixels that have the 
same sky positions to further remove fainter cosmic rays. The 
images we re calibrated by obse rving various UKIRT Faint 
Standards (lHawarden et alj|2001l) at least every half night on 
each detector. Data taken under nonphotometric conditions 
and poorly calibrated archive data were recalibrated with pho- 
tometric data taken by our group. 

The final mosaics have weighted exposure times (relative 
to median sky brightness and transparency) of 13.2 and 23.7 

8 also see http://www.aoc.nrao.edu/~whwang/idl/SIMPLE/index.htm 
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TABLE 1 
NIR Photometry of GOODS 
850-5 



Wave Band 


Flux (nJy) 


MOIRCS / 


-0.8 ±13.8 


NICMOS F160W 


-2.75 ±4.90 


MOIRCS Ks 


9.1 ±13.9 



hours in the J and K s bands, respectively, at the location of 
GOODS 850-5. The image qualities are very good, with 
FWHMs of 0."46 at J and 0."42 at K s . The rms astrom- 
etry error between the MOIRCS source positions and the 
ACS/SuprimeCam reference catalog is 0"08. The reduced 
MOIRCS images in the region around GOODS 850-5 are 
shown in Figure Q] The 1 a sensitivity limits at J and K s are 
both 14 nJy. As with the F160W image, GOODS 850-5 is not 
detected in the J and K s images. 

2.3. NIR Photometry 

Although GOODS 850-5 is not detected in any of the above 
ultradeep NIR images, it is useful to measure its fluxes at the 
SMA position (J2000 = 12:36:33.45, +62:14:08.65, W07) to 
determine whether there is any low level flux recorded. It 
is also important to obtain realistic flux limits/errors in these 
bands for the use of photometric redshift fitting. We placed 
apertures with diameters that are 1.5 times the FWHMs of 
compact objects in the fields (0"46, 0."3, and 0."42 in J, 
F160W, and K s , respectively) at the SMA position to measure 
the fluxes of GOODS 850-5. For each image, the background 
was estimated in a 3" area around GOODS 850-5 after de- 
tected objects and GOODS 850-5 itself were masked. The 
measured / and F160W fluxes are negative, and the K s flux 
is 0.65 a. To estimate flux uncertainties, we first masked all 
detected objects in the images and then randomly placed the 
apertures in the neighborhood of GOODS 850-5. The dis- 
persions in the random aperture fluxes are considered as flux 
uncertainties in these three bands. This procedure is neces- 
sary especially for the two MOIRCS bands, where the images 
were not drizzle combined and therefore the noise is more 
correlated between pixels. Such flux errors also include the 
uncertainties in background subtraction and confusion noise 
from faint undetected sources in the field. We summarize the 
measurements in Table Q] 

2.4. Existing Data 

GOODS 850-5 is not detected in the GOODS-N HST ACS 
images but is clearly detected in the Spitzer IRAC 3.6-8 /xm 
and MIPS 24 /im images (GOODS Spitzer Legacy Program, 
M. Dickinson et al. 2008, in prep aration). It s IRAC and 
MIPS fluxes were first measured by iPope et alj (12006b . Be- 
cause GOODS 850-5 is blended with brighter nearby IRAC 
and MIPS sources, in W07 we remeasured the IRAC and 
MIPS fluxes with a PSF fitting method to better isolate its 
flux. In this work, we adopt the W07 results. (The flux errors 
in W07 and here include the uncertainties in the deblending 
processes.) GOOD S 850-5 is not detected in the MIPS 7 /im 
and 160 Lim bands dFraver et al.l2006l:lHuvnh et al.l2007l) . We 
adopt the 1 a limits of 2 mJy at 70 /im and 30 mJy at 160 /im. 

In the submillimeter, we adopt the SCUBA 850 /im jig- 
gle map flux of Wang et al. (2004) and the SMA 880 /mi 
flux of W07. GOODS 850-5 was detected by D08 with the 
IRAM Plateau de Bure Interferometer at 1.25 mm with a flux 
of 5.0 ± 1.0 mJy. In W07, we used the older Very Large 



TABLE 2 

Photometric Data of GOODS 850-5 



Wave Band Flux (/tjy) Reference 



ACS F435W -0.013 ±0.004 2 

ACS F606W -0.004 ±0.003 2 

ACS F775W 0.001 ±0.006 2 

ACS F850LP -0.009 ±0.009 2 

MOIRCS J -0.0008 ± 0.014 1 

NICMOS F160W -0.002.8 ± 0.0049 1 

MOIRCS K s 0.0091 ±0.014 1 

IRAC 3.6 pirn 1.14±0.14 2 

IRAC 4.5 /im 1.64±0.13 2 

IRAC 5.8 /im 2.33 ±0.24 2 

IRAC 8.0 /im 5.37 ±0.37 2 

MIPS 24 /im 46.3 ±9.2 2 

MIPS 70 /im < 2000 3 

MIPS 160 /im < 30000 4 

SCUBA 850 /im 12900 ±2100 5 

SMA 880 /tm 12000 ± 1400 2 

IRAM 1.25 mm 5000 ±1000 6 

VLA20cm 34.4 ±4.2 6 



REF ERENCES. — (1) th is wor k; (2) IWang et all 

(20071 W07); (3) Frav er et alj f2006h : (4) IHuvnhetalJ 

120071); (5) IWang. Cowie. & Bargeri 120041) ; (6) 
IDannerbauer. Walter. & Morrisori <2008L D08V 

Arr ay (VLA) 20 cm image of iBiggs & Ivisonl (I2006L also 
see lRichardsll2000l) . in which GOODS 850-5 has a flux of 
18.7 ± 8.0 /iJy 9 . G. Morrison et al. (in preparation) obtained 
a deeper VLA 20 cm image of GOODS-N. GOODS 850-5 
was detected with a highly significant flux of 34.4 ±4.2 in the 
new VLA image (D08). Here we adopt the latest D08 values 
at 1 .25 mm and 20 cm. The photometry used in this work is 
summarized in Table Figure |2] show the multicolor images 
of GOODS 850-5. Comparing to the figure shown in W07, 
the main difference is the new Subaru and HST NIR images 
and the new VLA 20 cm image. 

3. SED AND REDSHIFT OF GOODS 850-5 

3.1. Non-detections in J, F160W, and K s 

The results of our new NIR imaging are surprising and not 
particularly easy to understand. GOODS 850-5 is clearly de- 
tected in the IRAC 3.6 /im band at 8 a. Despite the great 
sensitivities at F160W and K s , which are respectively 29 and 
10 times higher than those at 3.6 /im, GOODS 850-5 is not 
detected. This reveals an extraordinarily red SED in the NIR. 
The 1 a K s flux upper limit (23 nJy) implies a Si^^m/SKs flux 
ratio of > 50 or a spectral slope of a < -7.8. For comparison, 
objects with R-K > 5, corresponding to a < -2.5, are called 
"extremely red objects" (EROs, e.g., lMcCarfhyll2004l) . The 
spectral slope of GOODS 850-5 is more than three times that 
of EROs, and the K s -3.6 /im color (in AB scale) of GOODS 
850-5 is > 2 higher than the I — J color of Hu-Ridgway 10, 
which is a prototype dusty starburst ERO (Hu & Ridgway 
1994). High redshift UV emitting objects are commonly se- 
lected with red colors of > 2 (in the AB system) between 
two adjace nt optical filter bands (the "Lym an-break" tech- 
nique, e.g., ISteidel et alJll999t iFan et al]l2001h . correspond- 
ing to a < -10. The flux ratio of GOODS 850-5 between K s 
and 3.6 /im is ^ 8 times that of a minimal Lyman break. (A 
Lyman break between K s and 3.6 /im would imply z > 17.) 

9 This flux was measured at the location of the SMA position, which is 
~ 0"3 away from the best-fit VLA position. G. Morrison (2008, personal 
communication) and R. Ivison (2008, personal communication) found best- 
fit VLA_fluxej>_of24 ; _l± 5.8 and 23.8 ± 5.9 /tjy, respectively, from the image 
of lBiggs & Ivlsof] <2006|) . 
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FIG. 3 . — NIR colors of GOODS 850-5 and galaxies between z = an d 10: (a) colors between the F160W and the 3.6 /im bands; (b) colors between the K s and 
3.6 [im bands. Solid curves are colors of E, Sbc, and Im galaxy types in Coleman, Wu, & Weedman 1 1980), without reddening. Dotted curves are the same types 
of galaxies with Ay = 2.0 (a) and Ay =5.0 (b). We adopt the extinction law in lCalzetti et all 120001) . Horizontal dashed lines are the lower limits for the colors of 
GOODS 850-5, derived from the observed Icr upper limits of its F160W and K s fluxes. 

The extremely red color of GOODS 850-5 sets a strong con- 
straint on its redshift. It can be seen in Figure [3]that it is very 
difficult to explain the observed colors between 1.6 and 3.6 
/im with galaxies at z < 5 even with a considerable amount 
of reddening. Simply based on the two colors in Figure [3] 
z ~ 6 appears to be the most likely redshift for GOODS 850- 
5. However, we can use all the filter bands simultaneously to 
obtain better redshift estimates with the photometric redshift 
technique. 



3.2. NIR Photometric Redshift 

In W07 we derived photometric redshifts for GOODS 
850-5 and found that galaxies at z ~ 6 provided the best 
fits. Here we used our new photometric data in the NIR 
to improve our photometric redshift estimates. We used 
all the Spitzer IRAC, HST ACS and NICMOS, and Subaru 
MOIR CS data and the latest version o f the Hyperz pack- 
age dBolzonella. Miralles. & Pellol l2000h 10 for this calcula- 
tion. In each band, we added a 5% of zeropoint error 
quadratically. We used zeros to replace negative fluxes in 
the four ACS bands and the J and F160W bands, since 
negative fluxes are not physical. (The F775W flux is con- 
sistent with zero.) We used the non-zero K s flux. The 
photometric redshifts were derived independently with two 
SED template sets. The first is the lat e st stel lar population 
synthesis model of iBruzual & Charlotl (120031 BC03). The 
second contains the empiri c ally o bserved galaxy spectra of 
IColeman. Wu. & Weedmanl d 19801) f r om E to Im types and 
the starburst spectra of lKinnevetal.1 (119961) . Both SED sets 
are widely used by the optical extragalactic community for 
photometric redshift estimates. The combination of the latest 
Hyperz with the BC03 model has the nice feature of estimat- 
ing ages and stellar masses of galaxies. Hyperz also limits the 
age of the galaxies not to exceed the age of the universe. 

We adopted the standard extinction law of ICalzetti et al.l 
(2000). There is a subtlety in the maximum extinction to be 
allowed. Explaining the extremely red color of GOODS 850- 
5 between F160W//TJ and 3.6 /im requires a strong Balmer 
break (rest wavelength 4000 A) at a high redshift of z > 
5 (or a Lyman break at z > 16). However, if we allow 
very large extinctions, this red color might be reproduced 

10 also see http://www.ast.obs-mip.fr/users/roser/hyperz/ 



with a highly reddene d low-redshift galaxy. For example, 
Moba sher et al.1 (120051) suggest that the very red z — J color 
of the object HUDF-JD2 may come from a Lyman break at 
z ~ 6.5. On the othe r hand, the mid-IR (MIR) observations 
of IChary et all (120071) suggest that HUDF-JD2 is at z ~ 2 and 
the red z—J color comes from an extinction of Ay ~ 4.9 (als o 
see lFontana et al.l2"00 6; Dunlo pTCirasuolo. & McLurel 2007). 
Spectroscopical ly confirmed SMGs have a typical exti nction 
of Ay ~ 1-3 dSmailetal] 12004 ISwinbank etal] 120031) . but 
the spectroscopic sample may be biased toward less reddened 
systems. Photometric redshift fitting to SMGs without spec- 
troscopy generally gives Ay ~ 0-5. From these examples, 
Ay < 6 seems to be a reasonable limit for high-redshift SMGs. 
(For heavier extinctions, the shape of the extinction curve be- 
comes a much more important uncertainty than the Ay value 
itself.) We adopted this limit for our photometric redshift 
analyses. Later we will show that, unlike in the case of 
HUDF-JD2, a z < 3 redshift with a large extinction is ruled 
out for GOODS 850-5 by its radio and FIR SED. 

Figure|4]shows the results of our photometric redshift anal- 
yses. Both SED template sets give best fits at redshifts be- 
tween 6 and 7. The BC03 best-fit is a galaxy at z = 6.9 that 
formed in a single burst of star formation and has now aged 
0.7 Gy with an extinction of Ay = 1.8. The empirical SED 
set provides a best fit at z = 6.4, with an elliptical type and 
Ay = 0.5. The BC03 set systematically gives higher extinc- 
tions at all redshifts, likely because of the limit in the galaxy 
age, i.e., an unreddened galaxy that is sufficiently red to fit 
the observed SED may be too old for the age of the universe. 
Nevertheless, the xl distributions from the two SED sets are 
consistent with each other. In our subsequent analyses, we 
adopt the BC03 result as it includes an age limit for the galaxy 
and provides physical quantities such as the age and stellar 
mass. The 68% confidence range (1 a for a Gaussian distri- 
bution) from the integration of the BC03 probability function 
is z = 6.0-7.4. We note that redshifts between 4 and 5 are not 
entirely ruled out, although the fit is poorer here. On the other 
hand, z < 3 is safely ruled out by the strong limits from the 
NICMOS and MOIRCS non-detections. 

Comparing to the photometric redshift analyses in W07, the 
probability distribution now becomes narrower, which is an 
important improvement of this work. On the other hand, the 
minimum xt increases, indicating that it is generally difficult 
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FIG. 4. — NIR photometric redshift results. The top panel shows the ob- 
served SED between 0.4 and 8 fim, and the best-fit SEDs from the two SED 
template sets. Filled symbols are detections and the error bars are smaller 
than the symbols. Open symbols with arrows are 1 a upper limits. The open 
symbol with an error bar is the 0.65 cr S" s -band measurement. Note the ver- 
tical scale of the top panel. The bottom panel shows the distributions of the 
xt. Ay, and probability of the fits. 

to explain the red color between 1.6 and 3.6 /im. We tested 
photometric redshifts at z > 10 and obtained nearly perfect 
fits at z > 17 (i.e., a Lyman break between K s and 3.6 /im). 
However, this high redshift cannot explain the observed SED 
in the FIR and radio. We do not consider z > 10 in this paper. 

3.3. Millimetric Redshift 

We used the radio and FIR portion of the SED for another 
photometric redshift estimate (the millimetric redshift esti- 
mate). Because there is not a clear correlation between the 
radio/dust SEDs and the stellar SEDs of galaxies, this mil- 
limetric redshift estimate was carried out independently of 
the above NIR photometric redshift. The earliest version of 
such millimetric redshift estimates was carried out using the 
spectral ind ex between two wav e bands in the radio and sub- 
millimeter dCarilri & Yunl 119991; iBarger. Cowie. & Richards! 
120001: 1 Yun&Carillil 120021) based on the well known radio- 
FIR correlation in the local universe. This method has larger 
errors caused by the uncertainty in the dust temperature. Ad- 
vanced versions utilize full SED fitting in the radio and FIR 
as the amount of available data increases. 

Here we used all the data listed in Table [2] between 24 /im 
and 20 cm, including the two non-detections at 70 and 160 
/im. These two non-detections play a key role in ruling out 
low redshifts of z < 3 but do not provide strong constraints at 
z > 4. We used two sets of radio and FIR SED templates. 
The first set contains model SEDs of Arp 220 (ultralumi- 
nous starburst with cooler dust emission), NGC 6090 (lumi- 
nous starburst), M 82 (low luminosity starburst with warm 
dust) and M 100 (normal spiral) adopted from ISilva et aT] 
( 1998), and Mrk 231 (ultraluminous dusty active galactic nu- 
cleus, AGN, with warm dust) derived from the photometry 
in NED 11 The second SEP set includes all of the 105 SED 
models in lCharv & Elbazl (120011 CE01). 

11 The NASA/IPAC Extragalactic Database (NED) is operated by the Jet 
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FIG. 5. — Radio and FIR photometric redshift results. The top panel shows 
the observed SED of GOODS 850-5 and the best-fit SEDs. The bottom panel 
shows the reduced \ 2 of the fits vs. redshift. In the bottom panel, the black 
curve only shows the best-fit among all 105 model SEDs in CE01 at each 
redshift. 

In Figure [5] we present the best-fit SEDs and the redshift 
distribution of xt for eacn SED template. A low redshift of 
z < 3 is ruled out and this agrees with the NIR photometric 
redshift. The best-fit with the Silva et al. models comes from 
Arp 220 at z = 3.9. This is perhaps not too surprising since 
Arp 220 has been commonly considered as the local analogue 
to high-redshift SMGs. The next best fit comes from M 82. 
It is at z = 6.2 and the fit is slightly poorer than that with Arp 
220. The fit with the spiral template M 100 has xl > 10 ev- 
erywhere and is not shown in the figure. The best-fit AGN 
template has a redshift of 9.1, which was discussed in W07. 
The best-fit of the CE01 templates is at z = 4.5 and has a xi 
(~ 0) lower than that from the Silva et al. models. This un- 
reasonably low xt is perhaps a result of the very wide range 
of models (105 of them). Generally speaking, the CE01 mod- 
els provide good fits in the z ~ 3.5-5.5 range. D08 used both 
the two-waveband spectral index method and the SED fitting 
method (also with the CE01 templates, but only with the 24 
/im, 850 /im, 1 .25 mm, and 20 cm bands) and obtained z ~ 3 .7 
and z ~ 3.3 respectively. Our photometric redshift results are 
slightly higher than those in D07, likely due to the introduc- 
tion of the two upper limits at 70 and 160 /im, which strongly 
disfavor z < 3.5. Nevertheless, the results here and those in 
D08 are still broadly consistent. 

3.4. z~4orz>6? 

We are now faced with two different photometric redshift 
results. The NIR photometric redshift suggests z ~ 6-7, while 
the millimetric redshift suggests z ~ 4. We favor the NIR red- 
shift of z > 6 because we believe that our knowledge about 
the stellar photospheric emission of galaxies is more robust 
than that about the dust and radio emission. This is not 

Propulsion Laboratory, California Institute of Technology, under contract 
with the National Aeronautics and Space Administration. 
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only because the stellar spectral synthesis models (such as the 
Bruzual & Chariot models) and the empirical spectra (such as 
the Coleman, Wu, & Weedman templates) have been widely 
adopted and tested in numerous extragalactic studies from 
z ~ to z > 7, but also because of the various uncertainties 
in the radio/IR photometric redshift. Below we discuss each 
of the uncertainties. 

First, the two-waveband spectral index method suffers from 
the uncertainty in dust temperature. Second, the stellar con- 
tamination in the 24 /im band is quite uncertain at z > 4. All 
NIR photometric redshift models for GOODS 850-5 at z > 4 
predict observed 24 /im stellar photospheric fluxes that are 
much weaker than those in the Silva et al. and CE01 models. 
If we decrease the stellar emission in the Silva et al. models 
to match the observations, all fits become better than those in 
Figure [5] In particular, the M 82 fit at z ~ 6 slightly shifts to 
z ~ 5 .5 and becomes better than the Arp 220 fit at z ~ 4. Third, 
and most importantly, there is a fundamental problem in the 
millimetric redshift, which is the uncertainty in the radio-FIR 
correlation. 

The radio-F IR correlation is a tight correlation observed for 
local galaxies dCondonl 19921) . but the detailed physical mech- 
anism of this correlation is still unclear. Several groups have 
attempted to measure this correlation on different samples of 
high-redshift galaxie s with various methods, and the results 
remain controversial ( Aroleton et al. 2004; Boyle et alJl2007h 
IVlahakis. Eales. & Dunndl2007h . If high-redshift SMGs are 
brighter in the radio, as suggested by Vlahakis et al., then 
the redshift of GOODS 850-5 should be higher than that in- 
ferred from the local radio-FIR correlation. In addition to the 
uncertainties in observational determinations of the correla- 
tion, there are also physical reasons to suspect that the radio- 
FIR correlation at high redshift would be different than that 
in the local universe. For example, the higher energy den- 
sity of the cosmic microwave background radiation at high 
redshift leads to a str onger inverse -Compton cooling for the 
cosmic-ray electrons (ICondonlll992h . which suppresses the 
radio emission. If this is the dominant effect at high redshift, 
the correct redshift of GOODS 850-5 should be lower than 
that inferred from the local correlation. On the other hand, 
galaxies tend to have smaller sizes at high redshift. This may 
enhance the interstellar magnetic field and the synchrotron ra- 
diation, but the total radio emission can be further compli- 
cated by ionization and bremsstrahlung loss (iThompson et alJ 
l2006h . 

Unfortunately, at least in the case of GOODS 850-5, the 
result of the millimetric redshift highly relies on the validity 
of the radio-FIR correlation. To show this, we repeated the 
same photometric redshift fitting but did not use the 20 cm 
data point. The resultant \l distributions are shown in Fig- 
ure [6] In this experiment, all the fits have \l around 1.0 (the 
expectation value) between z ~ 3 and z > 10. This shows 
that the redshift of GOODS 850-5 highly depends on its dust 
temperature when the radio-FIR correlation is not assumed. 
While it is established that many high-redshift SMGs have 
cool dust SEDs similar to that of Arp 220, a small number 
of SM Gs have warmer, M 82-like SEDs (e.g JClements etalJ 
2008, also see Figure 2 in W07). Because of this, we cannot 
rule out higher redshifts of z > 6 and an M 82-like SED for 
GOODS 850-5. 

3.5. Summary on the Redshift 

A low redshift of z < 3 is ruled out by both the NIR and 
millimetric redshifts for GOODS 850-5. The NIR photomet- 
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FIG. 6. — FIR-only photometric redshift results. Description the curves 
is as in Figure [5] The plot legend (except for CE01) are ordered by dust 
temperature: Arp 220 has the coolest dust and Mrk 231 has the warmest dust. 
Note that the vertical scale is different than that in Figure|5] This plot shows 
that the observed FIR SED of GOODS 850-5 can be reasonably fitted with 
various templates from z ~ 4 to 10. 

ric redshift suggests z ~ 6-7, but the millimetric redshift sug- 
gests z ~ 4. The later strongly relies on the validity of the local 
radio-FIR correlation. If we do not assume the local radio- 
FIR correlation, redshifts between 3 and 10 seem equally pos- 
sible for the observed FIR SED of GOODS 850-5. We con- 
sider the z > 6 NIR photometric redshift as a more likely one, 
but z ~ 4 is still a possibility for GOODS 850-5. 

4. PROPERTIES OF GOODS 850-5 

4.1. Active Galactic Nucleus? 

The possibility of the existence of an AGN was discussed 
by W07. T he lack of an X-ray cou nterpart in the 2 Ms Chan- 
dra image (lAlexander et al.ll2003l) rules out a Compton-thin 
AGN at z < 3. An AGN with an X-ray luminosity ~ 10 42 erg 
s" 1 at z > 3 cannot be detected by Chandra and therefore is 
still possible. If the radio and FIR emission of GOODS 850-5 
is entirely powered by a Compton-t hick dusty AGN like Mrk 
231, its redshift would be z ~ 9 (§[33] W07), which is less 
likely. Furthermore, our new observation in the NIR shows 
a strong spectral break around 2-3 /im, which is an impor- 
tant feature of old stellar populations. On the other hand, 
QSOs typically have power-law SEDs across a wide wave- 
length range from the UV to the NIR and do not show strong 
spectral breaks like this. We conclude that there is no evidence 
for an AGN in any of the observations. 

4.2. Stellar Population 

The use of the BC03 models and the latest Hyperz allows us 
to estimate the stellar masses and ages. The SED of GOODS 
850-5 is nearly a power law in the IRAC bands but shows a 
clear spectral break between 1.6 and 3.6 ^m. This can hardly 
be explained by a normal extinction curve. A spectral break in 
the unreddened stellar SED is required, and the most natural 
spectral break is the 4000 A Balmer break, which is a signa- 
ture of older stars. Hyperz agreed with this interpretation and 
all the best-fit models between z = 4 and 10 have nearly the 
age of the universe. For example, at the best fit of z = 6.9, the 
age of the universe is 0.79 Gyr and the best-fit model has an 
age of 0.7 Gyr . The minimum age with good fits (xl proba- 
bility > 0.5 x the best-fit one) is 0.5 Gyr, corresponding to a 
formation redshift of 10-14. Because there are no detections 
in the K s and F160W bands, the actual strength of the Balmer 
break in GOODS 850-5 is unknown. Therefore we do not 
think we can constrain its age with sufficient accuracy. Nev- 
ertheless, we can conclude that it requires a well established 
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FIG. 7. — Stellar mass vs. redshift, derived from Hyperz with the BC03 
models. Solid curve is the best-fit result. Dotted curves are the maximum and 
minimum stellar masses when the \v i s larger than the minimum xi by 1.0 
(see Figure|4]for the distribution of minimum xJX showing how the stellar 
mass changes when the fit is perturbed around the best-fit. At z < 4, the fits 
become very poor although the best-fit mass still follows the trend at z > 4. At 
z ~ 6.9, the minimum and maximum range is approximately ±40% around 
the best fit. There is another factor of 2 uncertainty caused by the uncertainty 
in the stellar population model (see text). 

stellar population with an age that is significantly large com- 
pared to the cosmic time to explain the observed SED at < 10 
/im. 

Since the observed IRAC fluxes are likely dominated by old 
stars with little AGN contributions, they can be used to make 
a measurement of the stellar masses. Figure [7] shows the Hy- 
perz/BC03 best-fit masses as a function of redshift, and the 
uncertainty range. At z = 4 and 6.0, the best-fit masses are 3 
and 5 x 10 11 M Q , respectively. The range allowed by the pho- 
tometric redshift fitting is quite narrow at z > 4. We therefore 
believe that the true uncertainty in the mass estimate is more 
likely from the stellar population model. In BC03, we can see 
that the mass-to-light (M/L) ratio is a strong function of ini- 
tial mass function (IMF) and a w eaker function o f metallicity. 
The adopted IMF in BC03 is the lChabrieri(l2003l) IM F, which 
is more top-heavy than the standard ISalpeterl d 1955b IMF. If 
the IMF is Salpeter, the stellar mass of GOODS 850-5 would 
be approximately 2 times larger than the above values. An- 
other important factor is the likely existence of the thermally 
pulsing asymptotic giant branch (TP-AGB, Maraston 2005) 
phase, which is not included in the BC03 models. Such a 
phase decreases the M /L ratio for a yo ung stellar population 
in the rest-frame NIR. iMarastor] ([2006) found that the inclu- 
sion of the TP-AGB stars decreases the mass by ~ 60% on 
average, compared to the BC03 models. This would decrease 
the stellar mass of GOODS850-5 to 3 x 10 11 M Q at z - 6, but 
this is still a very large mass. 

4.3. IR Luminosity and Star Formation Rate 

The data at > 24 /im allow us to robustly determine 
the total IR luminosity (integrated from 8-1000 /im, e.g., 
ISanders &~ Mirabel 1996). Figure [8] shows the IR luminosity 
of the best-fit ISilva et alJdl998b models (FIR only, §[34] and 
Figure[6]) as functions of redshift and SED type. We consider 
two redshifts here: z ~ 4.0, as suggested by the millimetric 
redshift estimate, and z ~ 6.9, as suggested by the NIR photo- 
metric redshift. At z ~ 4 the best-fit model is Arp 220 and the 
best-fit IR luminosity is 1.4 x 10 13 L Q . At z ~ 6.9 the best-fit 
model is M 82 and the luminosity is 2.6 x 10 13 Lq. These 
are all comparable to the result directly inferred from the sub- 



er, & Cowie 



1 1 1 1 1 


- ■** — 


" 






- 












Arp220 




N6090 




M82 


/ I 


Mrk231 


f; 




/•' 




/•' 




'.■ 

' ' ■ ■ 





2 3 4 5 6 7 8 
Redshift 

FIG. 8. — Total IR luminosity for the best-fit models vs. redshift. The IR 
luminosity is derived from SED fitting to data between 24 fim and 1 .25 mm. 
The result from the M 100 SED is not shown since M 100 does not fit the 
observations (xj > 10 for all redshifts). 

millim eter flux using the standard formulas (e.g., Blain et al. 
120021 L IR = 1.9 x 10 12 Sg 50f , m Lo/mJy, which is 2.5 x 10 13 L 
for GOODS 850-5). 

The star formation rate for GOODS 850-5 can be estimated 
with the relation M = 1.7 x 10" 10 L IR /L Q (e.g., iKennicutj 
119981) . The results are 2400 M yr" 1 for z - 4 and 4400 M 
yr" 1 for z ~ 6.9. This assumes the standard Salpeter IMF. 

It is interesting to compare the above star formation rates 
with the non-detections in the NIR. Using t he UV luminos - 
ity versus star formation rate conversion in iKennicuttl (1 19981) . 
M = 1.4 x 10" 28 Luv/(erg s" 1 Hz" 1 ), and a minimum IR star 
formation rate of 2000 M Q yr" 1 , we found an unattenuated 
UV (1500-2800 A) luminosity of 1.4 x 10 31 erg s" 1 Hz" 1 . At 
z = 4-6.9, this corresponds to approximately 40-20 /iJy at 1.6 
to 2.1 /im, assuming a young burst SED from BC03. The ob- 
served F160W and K s flux upper limits then imply the extinc- 
tions for the starburst component to be at least Ay ~ 6.7-4.6. 
Figure [9] shows the SEDs of such hidden starbursts at z = 4 
and 6.9. The lower limits for the extinctions of the starburst 
component are significantly higher than those derived from 
the optical/NIR SED. We will come back to this in § 1431 

4.4. Nature of the 24 fim Emission 

All the above discussions on the existence of an AGN, the 
stellar mass, and the IR luminosity depend (to various de- 
grees) on whether the observed 24 /im emission is dominated 
by emission from stellar photosphere or from dust (see the 
e xtensive discussion on the 24 /im emission from HUDF-JD2 
in lMobasher et alj|2005l) . Roughly speaking, the observed 24 
/im flux of GOODS 850-5 is unlikely to be dominated by stel- 
lar emission, as this would require Ay > 8 at z ~ 6 to explain 
the observed red 24-to-8.0 /im color without dust emission. 
If we include the 24 /im band in the optical/NIR photometric 
redshift analyses, the fit becomes poor at all redshifts. 

We can quantify the strength of stellar radiation at 24 /im 
using photometric redshift fits without including the 24 /im 
data. The best-fit BC03 models in § 13.21 imply 24 /im stel- 
lar emission of 4.6 and 2.9 /iJy at z ~ 4 and 6.9, respectively, 
corresponding to 10% and 6% of the observed 24 /im flux. 
These seem unusually small, especially at z ~ 6.9, where 24 
/im corresponds to res t-frame 3 /im. We can compare these 
with local templates in lSilva et alJ (119981) to see if they make 
any sense. Table[3]compares the ratios of Lm/M+ and the ra- 
tios of dust-to-stellar radiation at 24 /im for Arp 220, NGC 
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FIG. 9. — NIR SEDs of ongoing starbursts suggested by the IR luminosity 
of GOODS 850-5 at z, = 4.0 and 6.9. Two curves are shown for each redshift: 
one with no extinction, derived from a BC03 instantaneous burst model with 
an age of 5 Myr and a star formation rate of 2000 Mq yr ; the other with an 
extinction that is required to minimally hide the burst in the F160W and K s 
bands such that the hidden burst contributes < 50% to the observed 1 a upper 
limits. The squares are the observed SED of GOODS 850-5 at > 1 jum. 



TABLE 3 

Dust and Stellar Contributions at 24 fiM 



Z Galaxy 




c-dust / c* 
J 24 l°2A 




4 


Arp220 


14 


4.8 


0.34 


NGC 6090 


1.5 


1.6 


1.1 


M82 


0.2 


2.7 


13 


GOODS 850-5 


46 


15 


0.33 


6.9 


Arp 220 


14 


0.08 


0.0057 


NGC 6090 


1.5 


0.07 


0.047 


M82 


0.2 


0.08 


0.40 


GOODS 850-5 


52 


9.1 


0.18 


Note. — ^ir/M* 


is in solar 


unit. 


(S^>/S* 4 )' = 


(S% s '/S* M )/(L m /M+). 









6090, M 82, and GOODS 850-5, at z ~ 4 and z ~ 6.9. We 
found that once the dust-to-star ratios in the observed 24 /im 
emission is normalized by the Lir/M+ ratios, GOODS 850-5 
is not different from local galaxies. It is also interesting to 
note that at z ~ 4 and 6.9 the normalized dust-to-star ratios 
are similar to Arp 220 and M 82, respectively. This is consis- 
tent with what we found in the photometric redshift analyses: 
GOODS 850-5 is similar to Arp 220 if it is at z ~ 4 but more 
similar to M 82 if z > 6. 

There are two subtle issues in the above comparison. First, 
the underlying assumption for the normalization is that the 
M/L ratios are similar for the galaxies at 24 /im (rest-frame 
3.0^4.8 /im for the redshifts of interest). Given the long wave- 
lengths, this is acceptable. Second, the estimates of the IR 
luminosity for GOODS 850-5 slightly depends on how much 
of the observed 24 /im emission comes from the warmer dust 
components. However, from Figure[8]we can see that the un- 
certainty is well within a factor of 2 and the IR luminosity 
is mostly determined by the three measurements in the mil- 
limeter and submillimeter. Therefore the 24 /im measurement 
can still be considered as a semi-independent check for all the 
analyses in longer and shorter wavelengths. The fact that the 
extrapolations of the optic al/NIR photometric redshift analy- 
ses agree well with the millimetric one at 24 /im (in terms of 



dust-to-star ratio) does show an excellent self consistency. It 
is fair to conclude that the unusually large dust-to-star ratio in 
the observed 24 /im emission of GOODS 850-5 is simply a 
consequence of its extremely large IR luminosity. 

Two question arises once we establish that > 90% of the 
observed 24 /im emission comes from dust: (1) Does this 
affect our millimetric redshift analyses? (2) Does the dust 
emission extend to the IRAC bands and bias our stellar mass 
estimates? To examine (1), we decreased the stellar contribu- 
tion in the Silva et al. templates by factors of 2 to a few tens 
and repeated the photometric redshift fitting. We found that 
the fits improve slightly. The best-fit redshift with Arp 220 
does not change significantly and the best-fit redshift with M 
82 decreases from ~ 6-7 to ~ 5.5-6. Thus, the general con- 
clusions on the millimetric redshift and the IR luminosity are 
both fairly insensitive to the assumed stellar contribution to 
the 24 /im flux. 

The answer to (2) is negative as well. The dust spectral 
slope for starbursts is extremely steep at < 1.6 /im (IRAC 
bands for z > 4), and the blue end of the dust emission does 
not contaminate the IRAC fluxes. This is still true even if the 
observed 24 /im emission come s from a dusty AGN wi th dust 
much warmer than starbursts. Mobas her et alJ (120051) com- 
pared the observed 24 /tm emission from HUDF-JD2 with 
dusty AGN templates Mrk 231 and NGC 1068. We follow 
the procedure in Mobasher et al. and compare the SEDs of 
GOODS 850-5 and NGC 1068. We found that the dust emis- 
sion only contributes to the IRAC fluxes at z < 4. At z > 6, 
similar to what Mobasher et al. concluded for HUDF-JD2 at 
z = 6.4, the dusty AGN does not bias the IRAC stellar mass 
estimates. We do not consider Mrk 231 here as its SED in the 
literature is contaminated by stellar emission at rest-frame < 3 
/tm (as pointed out by Mobasher et al.) and hence it does not 
provide a definite answer. Moreover, as shown in our photo- 
metric redshift analyses, to explain the large millimeter fluxes 
of GOODS 850-5 with a Mrk 231 -like warm dust, it would 
require unlikely high redshifts of z ~ 8-10. In short, we con- 
clude that the dust emission at 24 /im is unlikely to bias our 
IRAC-based stellar mass estimates regardless of the existence 
of a dusty AGN. 

4.5. Coexistence of Old Stars and a Starburst 

The properties of GOODS 850-5 are puzzling. There ap- 
pear to be two inconsistencies that are related to each other. 
(1) While the IR luminosity implies an intensive ongoing star- 
burst and the formation of a young galaxy, the NIR SED 
suggests that most of the observed stellar radiation comes 
from old stars without any detectable rest-frame UV radia- 
tion. (2) Our interpretation of the NIR SED suggests Ay < 2, 
but Ay > 4.6 is required to hide all the young stars. 

One conceivable way to solve this is to see if Ay > 4.6 with 
young stellar populations could provide reasonably good fits 
in the photometric redshift. If yes, we may argue that the en- 
tire young galaxy is behind an Ay > 4.6 dust screen and that 
the observed IRAC fluxes indeed mostly come from young 
stars. To investigate this possibility, we limited the photo- 
metric redshift fitting to the BC03 models with Ay = 4.6-10, 
z > 4, and ongoing star formation. At z ~ 4, fits with such 
highly extinguished stellar populations all have xt > 2.7, 
which is significantly poorer than the fits with Ay < 2. At 
z > 5, the x„ becomes larger than 5 and the models fail to fit 
the data. We conclude that Ay > 4.6 for a young galaxy is a 
less likely case, although not fully ruled out for z ~ 4. It is 
more likely that the observed IRAC fluxes mostly come from 
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a relatively old stellar population with a moderate extinction 
ofA v <2. 

Another possible and indeed sensible scenario to explain 
the above inconsistencies is that the star forming region (with 
Ay > 4.6) in GOODS 850-5 is different from the region 
that produces most of the observed IRAC fluxes (Ay < 2). 
GOODS 850-5 is not resolved by the SMA (~ 2."2 beam 
FWHM, S/N=8.6) and by the VLA (~ l."7, S/N=8.2). These 
imply an uncertainty of ~ 0."2 (beam FWHM divided by S/N) 
for its beam-convolved size and ~ 0"8 for the upper limit of 
its intrinsic size. The upper limit corresponds to ~ 4.5 kpc 
at z = 6 .9 and is consisten t with the CO sizes of low-redshift 
SMGs (iTacconi et al.ll2006l) . On the other hand, the IRAC 
fluxes suggest a massive stellar population of ~ 10 1L5 M©, 
which may be spatially offset from the starburst region and 
may have less extinction. 

An interesting object for comparison is the prototypica l 
dusty starburst ERO, Hu-Ridgway 10. IStern et al.l (2006) 
found that its 10 /im silicate feature implies Ay ~ 11 but its 
optical/NIR SED implies a much smaller Ay ~ 2.4. They 
ascribed this to a heavily obscured starburst in Hu-Ridgway 
10, similar to our hypothesis for GOODS 850-5. Unfortu- 
nately, the current data of GOODS 850-5 do not allow us to 
use the silicate feature for an extinction measurement. An- 
other possibility to test the above two-component hypothesis 
for GOODS 850-5 is to see if we can fit the observed NIR 
SED with two components. We found that the photometric 
redshift fitting becomes poorer at all redshifts if we include 
hidden bursts shown in Figure [9] The increase in xt i s fr° m 
~ 0.2 at z = 4 to ~ 0.5 at z = 6.9. However, we note that the 
actual extinction for the starburst component can be higher 
than the minimally hidden burst shown in Figure [9] Once we 
increase the extinction for the burst component by ~ 2-3 (still 
less than the extinction in Hu-Ridgway 10), the photometric 
redshifts reduce to the ones shown in § !3.2l at all redshifts. In 
other words, the current data are insufficient for further testing 
the two-component hypothesis. 

We conclude that a massive galaxy with an old stellar 
population but with some dust screening hosting a much 
more dusty nuclear starburst can explain the observations of 
GOODS 850-5. 

5. DISCUSSION 

For the stellar masses and star formation rates estimated at 
Z ~ 4 and 6.9, the mass build-up times are both ^110 Myr, 
respectively corresponding to 7% and 14% of their cosmic 
times. It is clearly more difficult to produce the stellar mass 
for the z ~ 6.9 case, since there is less time available for the 
formation of the galaxy. Given that all the stars we see are 
old, it requires an intensive but short burst of > 10 3 M© yr" 1 at 
Z > 10. This also suggests that GOODS 850-5 has undergone 
two distinct bursts of star formation: one produced the old 
stars observed in the IRAC bands, and the other is producing 
the observed dust emission in the MIPS and (sub)millimeter 
bands. 

The estimated mas s of M+ ~ 10 11,5 Mp) is similar to t ypical 
SMGs at z < 4 (e.g. JSmail et al.ll2004t lDvelFaIll2008l) but is 
more massive than other galaxies observed at z > 6. The most 
massive Lyman-break (UV luminous) selected galaxies at 
z ~ 6-7 have stellar masses of ~ 10 10 — 10 10 5 M© (lLabbe et all 
120061 lYanetal.1 120061: lEvles et alJl2007h . GOODS 850-5 is 
10 times more massive than these optically selected galaxies, 
suggesting that the very massive galaxies at these redshifts 
may be still in dusty starbursting phases and therefore may be 



missed by deep optical surveys. 

We can estimate the limit for the number of massive, dust- 
hidden galaxies at z > 6. The photometric redshift uncertainty 
range for GOODS 850-5 is z ~ 6.0-7.4 and our GOODS- 
N SCUBA survey area for sources brighter than GOODS 
850-5 is ~ 100 arcmin 2 . Therefore the comoving volume 
of our survey for sources similar to GOODS 850-5 is ~ 
3 x 10 5 Mpc 3 . This implies a number density of ~ 3 x 10~ 6 
Mpc" 3 for massive objects of M+ ~ 10 1L5 M©. It is unclear 
how a radio-faint submillimeter selection is related to stel- 
lar masses. However, the above density is probably a lower 
limit, since GOODS 850-5 is the first radio-faint SMG stud- 
ied in our SMA survey and the survey is not yet complete. 
This lower lim it is slightly lar ger than the maximum value 
suggested by lYan et al.l (2006, see their Figure 11) for op- 
tically selected galaxies at z ~ 6, and is significantly larger 
th an the values from the AC DM hydrodynamic simulations 
of iNagamine et all (120041) and lNight et all (12006). This em- 
phasizes the importance of dusty galaxies at high redshift, but 
we clearly need a much larger high-redshift SMG sample to 
reach cosmologically meaningful conclusions. 

It is interesting to ask whether current ACDM models can 
produce 10 11 5 M© galaxies at all at z > 6. In standard ACDM 
models, massive galaxies form in later cosmic times (aka. "hi- 
erarchical" galaxy for mation). On the oth er hand, evidence of 
"cosmic downsizing" (ICowie et al.lll996l) has been observed 
at many different wavelengths over a broad range of redshifts, 
and this has been considered by some as anti-hierarchical. 
The existence of ultralumino us SMGs at z > 2 once put a 
strain on the ACDM models dBaugh et al.ll2005l) . Moreover, 
the stacking detection of a population of faint SMGs at low 
redshifts of z ~ 1, which dominates the submillimeter EBL, 
further demonstrated the do wnsizing behavior in the SMG 
population from z <y 4 to z = (Wa ng. Cowie. & Bargedl2006t 
Seriea nTet al.l20 08). Recent semi-analytic and hydrodynamic 
models are able to reproduce galaxies with M+ > 10 10 Mp, up 
to z ~ 5-6 under the ACDM framework dBower et alJl2006t 
iNieht et al.ll2006h . 

At slightly h igher redshi f ts of z = 6-7, the ACDM N-body 
simulations of Luki c et ail (120071) provide a halo density of 
- 10" 5 (Mpc/h)" 3 for the mass range of 10 12 -10 13 M©. The 
density rapidly drops to 10 -9 — 10 _1 1 (Mpc/h)" 3 for the mass 
range of 10 13 -10 I4 M©. If we assume a matter-to-star mass ra- 
tio of > 10, the density of ~ 3 x 10" 6 Mpc" 3 that we derived 
for high-redshift SMGs seems to match the upper bound in 
the simulations. This suggests that it is plausible to find ap- 
proximately one z > 6 dark halo in our survey area that is mas- 
sive enough to host a hyperluminous SMG. However, whether 
the simulations can reproduce the observed intensive starburst 
and the large stellar mass within such a halo in a short cosmic 
time remains to be seen. 

Lastly, can objects like GOODS 850-5 play a role in the 
reioni zation of the unive rse? The latest 5 yr WMAP re- 
sult (jDunkley et al. 2008) implies a reionization redshift of 
11.3± 1.4 (for instantaneous reionization) and the massive 
old stellar population in GOODS 850-5 implies intensive star 
formation at z > 10- We us e the standard formulation in 
iMadau. Haardt. & Reesl (fl999b . p = 0.013 x f£ x [(1 +z)/6] 3 
M© yr" 1 Mpc" 3 , where p is the minimum star formation rate 
density required for reionization, / esc is the escaping fraction 
of ionizing photons, and a Salpeter IMF is assumed. With the 
number density estimated above, / esc ~ 0. 1, and assuming that 
the initial burst in GOODS 850-5 is as intensive as the current 
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one, we found that the formation of GOODS 850-5 at z > 10 
contributed < 10~ 2 to the ionizing photons that are required 
for reioinization. The fraction would be even smaller if the 
initial burst of GOODS 850-5 were dusty (i.e., a smaller /esc)- 
This is consistent with the current picture that the universe is 
ionized by a large amount of low luminosity objects. 

6. SUMMARY AND FINAL REMARKS 

Our new ultradeep NIR observations reveal many unusual 
properties of GOODS 850-5. It is detected by the SMA, 
IRAM, and VLA from the submillimeter to centimeter wave- 
lengths, and by Spitzer between 3.6 and 24 /zm, all with high 
significance. On the other hand, it is not detected in the J, 
F160W, and K s bands even with the nano-Jansky sensitivities. 
We analyzed the photometric redshifts of GOODS 850-5. The 
NIR photometric redshift suggests a high redshift of z ~ 6.9 
and rules out z < 3. The millimetric redshift also rules out 
z < 3 and suggests z ~ 4 if we assume the local radio-FIR 
correlation. Without this assumption, z = 4 to 10 is equally 
possible for the observed IR SED. We conclude that z > 6 is 
more likely for GOODS 850-5 but z ~ 4 remains a possibility. 

To explain the observed NIR SED of GOODS 850-5 and the 
IR luminosity requires an established stellar population that 
formed at a large look-back time coexisting with an intensive 
ongoing starburst that is completely invisible in the rest- frame 
UV. The old stars observed in the IRAC bands have a large 
mass of M* ~ 10 1L5 Mq and are ~ 10 times more massive 
than optically selected massive galaxies at z > 6. The current 
burst of star formation seems to be compact with a total IR 
luminosity of 1.4-2.6 x 10 13 Lq and a star formation rate of 
2400-4400 M Q yr" 1 . It is deeply enshrouded by an Ay > 4.6 
dust screen so its UV radiation is not detected by NICMOS 
at 1.6 /im and by Subaru at 2.1 /im. This large extinction 
required for the starburst component also makes it difficult to 



directly confirm its existence in the IRAC bands and the two- 
population hypothesis remains to be tested. 

The high redshift of GOODS 850-5, if confirmed, will have 
important implications for galaxy formation and evolution, as 
discussed in this paper and in W07. Its extreme faintness 
at < 2 prevents any optical and NIR spectroscopy with 
current space and ground-based instruments. A precise mea- 
surement of its redshift will most likely come from millime- 
ter "redshift machines" on large telescopes such as the 110 
m Green Bank Telescope. Spectroscopy in the MIR and FIR 
with the Herschel Space Observatory is another possibility. 
In general, the observations of GOODS 850-5 remind us that 
there is a class of unexpected objects that are not included in 
the traditional picture of radio identified SMGs. The SMA 
is likely to reveal more examples in the near future. Studies 
of cosmologically large samples of such high-redshift SMGs, 
however, will require combinations of next-generation instru- 
ments, such as the Expanded VLA, Atacama Large Millime- 
ter/Submillimeter Array, James Webb Space Telescope, and 
large bolometer arrays on single-dish millimeter telescopes. 
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